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Natural biological systems have evolved a diverse array of switches to realize their strategies for 
environmental response and development. Emerging applications of synthetic biology have begun 
to exploit such switches to achieve increasingly sophisticated designed behaviors. However, not all 
switch architectures allow facile design of the switching and memory properties. Furthermore, not 
all designs are built from components for which large families of variants exist, a requirement for 
building many orthogonal switch variants. Therefore, there is a critical need from genetic engineers 
for scalable strategies that yield custom bistable switches. Here, we use a sigma factor and its 
cognate anti-sigma factor to experimentally verify that ultrasensitivity from sequestration 
combined with positive feedback is sufficient to build a bistable switch. We show that sequestration 
allows us to predictably tune the switching boundaries, and we can easily tune our switch to 
function as a set-reset latch that can be toggled between two states by a pulse of inducer input. 
Molecular Systems Biology 8:620; published online 23 October 2012; doi:10.1038/msb.2012.52 
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Introduction 

Bistable switches represent the canonical solution for building 
memory units that also generate noise-tolerant sigmoidal 
responses to inputs (Ferrell, 2002; Burrill and Silver, 2010). 
Bistability in a dynamical system, such as a genetic switch, 
requires positive feedback (either directly or indirectly, as shown 
in Figure 1A) and any source of ultrasensitivity that generates a 
sigmoidal response (Ferrell, 1996; Angeli et al, 2004; Ray et al, 
2011). While Figure IB shows a number of commonly proposed 
sources of ultrasensitivity, currently all synthetic bistable systems 
rely on cooperativity (Gardner et al, 2000; Becskei et al, 2001; 
Atkinson et al, 2003; Kramer and Fussenegger, 2005; Ajo- 
Franklin et al, 2007), with possible exceptions considered in the 
Discussion section. In addition, all known natural bistable 
systems also utilize cooperativity, although often in combination 
with other sources of ultrasensitivity (Ozbudak et al, 2004; Acar 
et al, 2005; Trunnell et al, 2011). 

Despite the prevalence of cooperativity, recent work has 
shown sequestration's importance in natural systems (Kim and 
Ferrell, 2007; Buchler and Louis, 2008), and sequestration-based 
ultrasensitivity has been demonstrated synthetically in eukar- 
yotic systems (Bashor et al, 2008; Buchler and Cross, 2009; Lu 
et al, 2011; Lee and Maheshri, 2012). Curiously, no one has yet 
demonstrated in vivo that sequestration with positive feedback is 
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sufficient to build a bistable switch, but the design shown in 
Figure 1C should have beneficial design properties. As the 
minimal requirements are an activator to build the positive 
feedback loop and a tight-binding anti-activator, the design is 
generalizable to many potential parts. In addition, Figure ID 
illustrates that the response curve from sequestration can be 
predictably adjusted by varying the expression level of the anti- 
activator, and we demonstrate how this property can be utilized 
to predictably tune the switching boundaries. 



Results 

We chose to demonstrate our sequestration-based switch in 
Escherichia coli [Ec], using the sigW sigma factor and its 
cognate anti-sigma factor (rsiW) from Bacillus subtilis as the 
activator and anti-activator (see Materials and methods for 
cloning details and Figure 2A for schematic). To build the 
positive feedback loop, we placed the sigma factor under 
control of a SigW-dependent promoter. Expression of the anti- 
sigma factor was controlled by an arabinose-inducible promo- 
ter. Finally, system output was monitored by placing YFP 
expression under control of the SigW-responsive promoter. 

We call the cells ON when the positive feedback loop is 
active and YFP is expressed, and cells without YFP expression 
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Figure 1 Required elements for building a bistable switch. (A) Positive feedback can be indirect from mutual inhibition or direct from autoregulation. (B) Four commonly 
proposed sources of ultrasensitivity. Cooperativity (shaded in red) is the most prevalent source in known bistable systems, but we demonstrate that sequestration 
(shaded in green) can be used to build a bistable switch with better design features. (C) The sequestration-based switch requires an activator that forms a positive 
feedback loop and an anti-activator that sequesters activator molecules. (D) The response curve from sequestration can be shifted to the left and right by changing the 
expression level of the sequestering molecule. 
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Figure 2 Sequestration with positive feedback is sufficient for bistability in vivo. (A) Schematic of construct implemented in APA4309 E. coii cells using a sigma factor 
and anti-sigma factor from B. subtilis. Cerulean (CFP) is expressed constitutively in all cells to detect dead cells and for normalization. (B) Cells previously grown in full or 
no arabinose induction were inoculated into intermediate concentrations and grown for 12 h. The right arrowheads in maroon color show the mean fluorescence of cells 
previously grown in no arabinose (initially ON). The left arrowheads in orange color show the mean fluorescence of cells previously grown in full arabinose induction 
(initially OFF). The y axis is the mean YFP of the population after normalizing each cell by its CFP fluorescence. The shaded area shows the hysteretic region where the 
system retains a memory of its previous state. (C, D) Histograms at indicated arabinose concentrations show bimodal distributions at the switching boundaries. The x axis 
shows log 10 YFP fluorescence normalized by each cell's CFP fluorescence and y axis shows the normalized cell count. At least 1 000 cells were imaged at each induction 
point. Source data is available for this figure in the Supplementary Information. 
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are considered OFF. The defining feature of bistability is a 
hysteretic region where the steady-state response of the system 
depends on the previous state. Supplementary Figure 1 
presents the nullclines for this system and illustrates why we 
expect to observe hysteresis as anti-sigma factor expression is 
varied. To experimentally locate this region, we grew cells in 
the absence of arabinose and at full induction (0.1%) 
overnight and then stored both sets of aliquots at - 80°C. 
The cells from the aliquots were then inoculated into 
intermediate arabinose concentrations and allowed to grow 
for 24 doublings before being imaged by fluorescence 
microscopy. Figure 2B shows the memory region where the 
system responds differently depending on whether it was 
previously ON or OFF. Cells that were previously in full 
arabinose induction remain OFF in the memory region and 
then begin to turn ON at 10 _5 % arabinose (orange arrow- 
heads facing left). Cells that were previously grown in the 
absence of arabinose are fully ON in the memory region and 
begin to switch OFF at 8*10 ~ 4 % arabinose (maroon-colored 
arrowheads facing right). Constitutive CFP expression was 
used to normalize the YFP expression and eliminate dead cells 
from the analysis. Figure 2B also shows our model overlaid 
with the experimentally determined points (see Supplemen- 
tary information for modeling details). Figure 2C and D show 
that the distribution of cells in the ON and OFF states shifts at 
the switching boundaries. We found that ON cells grew 
~4.5% slower than OFF cells (Supplementary Table 1). Using 
a model recently published for stochastic switching between 
phenotypes with different growth rates (Nevozhay et al, 2012), 
we found that ON cells switched OFF in the memory region 
at a rate of 0.0165 ± 0.014 per hour while OFF cells switched 
ON at a much slower rate of 6.1*10 _4 ±7*10 -5 per hour 
(Supplementary Figure 2). We also verified that the sigma 
factor does not exhibit intrinsic cooperativity (Supplementary 
Figure 3) and that the state of the switch does not 
affect expression from a constitutive Sigma 70 promoter 
(Supplementary Figure 4) . To demonstrate that our design is 
not specific to SigW, we also demonstrated hysteresis in a 
bistable system based on SigE and its cognate anti-sigma 
factor ChrR from Rhodobacter sphaeroides (Supplementary 
Figure 5) . 

To tune the switching boundaries of our sequestration-based 
switch, we added another transcriptional unit with the sigma 
factor under control of a promoter inducible with anhydrote- 
tracycline (aTc) (schematic shown in Figure 3A). In this 
system, the two induced proteins compete with each other. In 
the limit of high expression of anti-sigma factor and no sigma 
factor induction, all the cells are turned OFF. In the other limit 
where the sigma factor is strongly induced and no anti-sigma 
factor is present, all the cells are forced ON. In between the two 
extremes, we expect a bistable region whose boundaries vary 
with the total concentrations of the proteins. The numerically 
computed bifurcation diagram for this system is plotted in 
Figure 3B. We predict that cells at higher levels of sigma factor 
(aTc) induction require higher levels of anti-sigma factor 
(arabinose) induction to cross the boundaries. We experimen- 
tally verified this relationship by preparing aliquots of master 
cultures that were forced to be OFF or ON. Initially, OFF cells 
were grown at the inducer concentrations indicated in 
Figure 3B with left arrowheads, and we found that increasing 
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aTc concentration results in switching ON at higher arabinose 
concentrations, shown in Figure 3C. Similarly, initially ON 
cells were grown at the concentrations indicated with right 
arrowheads, and we found that increasing aTc concentration 
requires higher levels of arabinose to switch OFF, shown in 
Figure 3D. 

From the bifurcation diagram, we also predicted that we 
could design our switch to function as a set-reset latch. The 
latch is a true memory element rather than a hysteretic switch. 
It holds the state of the last input: set or reset. Figure 3E 
illustrates that the two inputs to the latch are arabinose and aTc, 
and Figure 3F shows the logic table the set-reset latch 
implements. To demonstrate the latch with APA4310 cells, we 
sequentially performed the steps illustrated in the bifurcation 
diagram of Figure 3G. We first initialized the latch into the reset 
state by forcing the cells OFF with full induction of anti-sigma 
factor (step 1 of Figure 3G and response shown in Figure 3H). 
We then inoculated the cells with both inputs set low, so the 
cells are in the memory state and remain OFF (step 2) . While 
maintaining the same concentration of arabinose, we pulsed 
the cells into the ON state by growing them in aTc and observed 
that upon inoculating them into the memory state with both 
inputs set low, the cells remained ON (steps 3 and 4) . Notably, 
the cells at step 4 are at the same inducer concentrations as 
used at step 2, but the cells after step 4 are ON compared with 
being OFF after step 2. After 12 h in the memory state, we 
observed that <2% of the cells switched to the incorrect state. 
Finally, we could reset the latch by pulsing the cells with excess 
arabinose to force them OFF again (step 5). 



Discussion 

While bistability from sequestration has been shown in vitro 
(Kim et al, 2006) , we have implemented a previously proposed 
design (Francois and Hakim, 2004) and verified in vivo that 
bistability can be generated without relying on cooperative 
transcription factors. Recent papers have suggested other 
methods for building bistability in the absence of cooperativ- 
ity. Modeling of the MprA/MprB system in mycobacteria 
illustrated that sequestration of a sigma factor by its anti-sigma 
factor could have a critical role in allowing bistability in the 
natural system, and they propose experiments to demonstrate 
their predictions in vivo (Tiwari et al, 2010). Alternatively, 
nonlinearity in growth rate combined with positive feedback 
could result in bistability (Tan et al, 2009), but the in vivo 
synthetic system did not exhibit steady-state hysteresis on the 
population level, indicating that building memory units based 
on this design could be problematic. Just recently, two positive 
feedback loops without cooperativity or other sources of 
ultrasensitivity were theoretically shown to be sufficient to 
generate bistability (Palani and Sarkar, 2011). While this is an 
elegant topology as the requirements are minimal and the 
design facilitates tuning, the in vivo system contained basal 
ultrasensitivity in one of the feedback loops, so this design 
remains to be demonstrated using pathways without inherent 
ultrasensitivity. 

Biological devices should be evaluated in terms of operating 
performance and scalability. The key operational feature of our 
design is that it allows predictable tuning of the switching 
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Figure 3 Sequestration-based bistability allows predictable tuning of the switching boundaries and design of a set-reset latch. (A) Schematic of system implemented in 
APA4310 cells with extra sigma factor production controlled by aTc induction (B) Bifurcation diagram overlaid with arrowheads denoting inducer concentrations where 
the system was investigated. Cells initially OFF were induced at the concentrations indicated with left arrowheads, and cells initially ON were induced at the 
concentrations marked with right arrowheads. (C) Histograms show that cells initially OFF switch ON at higher arabinose levels with increasing aTc induction. Cells were 
grown at the indicated arabinose and aTc concentrations (corresponding to left arrowheads from panel B) for 12 h, and at least 1000 cells were imaged per histogram. 
The percentage indicates the proportion of cells considered OFF (sum of bars in shaded area). (D) Cells initially ON require more arabinose to switch OFF with increasing 
aTc. Arabinose and aTc concentrations correspond to right arrowheads from panel B. (E) Schematic of set-reset latch. (F) Logic table corresponding to set-reset latch. 
(G) Using the previously determined switching boundaries, a set-reset latch can quickly be designed. Path numbers correspond to observed points in next panel. (H) The 
system retains a memory of the previous SET or RESET pulse. Steps 2 and 4 were analyzed after allowing the system to grow for 12 h, and cells were grown in the Set 
and Reset stages for 9 h. Source data is available for this figure in the Supplementary Information. 



4 Molecular Systems Biology 2012 



© 2012 EMBO and Macmillan Publishers Limited 



Sequestration-based bistability 
D Chen and A Arkin 



boundaries by adjusting expression of the sequestering 
protein. Another advantage is that the simple topology 
based on positive autoregulation lets us describe the key 
characteristics with analytical expressions (Supplementary 
information). Thus, if we desire a certain dynamic range or 
width of hysteretic region for a future application, the 
analytical results guide us in choosing optimal parameter 
values. One caveat with our system, however, is that the width 
of the hysteretic region is coupled with the tuning of the 
switching boundaries (Supplementary Figure 6 illustrates how 
the nullclines are shifted with increasing extra sigma factor 
production). 

The scalability of any biological device is limited both by the 
number of orthogonal parts and by host interaction. We chose 
to use an extracytoplasmic function (ECF) sigma factor as the 
ECF sigma factors and non-ECF sigma factors recognize 
different consensus promoters via distinct mechanisms of 
promoter binding (Lane and Darst, 2006). There are 43 groups 
of ECF sigma factors that recognize different consensus 
promoters, and Ec contains two ECF sigma factors that belong 
to different groups than SigW's group (Staron et al, 2009). As 
35 of the ECF sigma factor groups have cognate anti-sigma 
factors that should be orthogonal, the strategy presented here 
should be scalable. We observe a small amount of crosstalk 
between SigW and SigE and their cognate promoters 
(Supplementary Figure 7), so using these systems in the same 
cell could require understanding the determinants of ortho- 
gonality. Host interaction could also be problematic as we see a 
mild effect on the growth rate when the sigma factor is 
expressed, but these could be resolved in the future by 
reducing the maximal rate of sigma factor production from the 
feedback loop or by identifying and fixing the parasitic 
interactions with the host. 

Memory modules will be key components of synthetic 
biology applications in the future, so genetic engineers will 
need a large palette of different types to choose from. While the 
genetic toggle switch based on mutually inhibitory transcrip- 
tion factors has been the standard design for a set-reset latch 
(Gardner et al, 2000), each instance requires two cooperative 
transcription factors while our design requires only one 
activator/ anti-activator pair per switch. Future work could 
further compare different properties of the competing designs; 
for example, the toggle switch was shown to be more robust to 
growth rate changes than a positive autoregulatory loop 
(Klumpp et al, 2009) . An in vivo demonstration for a T-latch 
was also recently published, which features a push-on and 
push-off behavior (Lou et al, 2010). 

Although we have focused on bistability, memory units can 
also be built on recombinases (Ham et al, 2008; Friedland et al, 
2009), which may be especially attractive for long-term 
storage, but these devices lack a hysteretic region that provides 
buffering and the scalability of these parts is still unproven. A 
recent design based on recombinases especially demonstrated 
the stability of the device over 120 doublings, although 
switching was only 85% complete after a RESET pulse 
(Bonnet et al, 2012). In contrast, our device is 99% complete 
in switching after either a SET or RESET pulse of 9 h, but we 
find that the half-life of the ON state in the memory region is 
only 24 h. Our current device may not be suitable for 
applications that require long-term stability, but this could 
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be improved in the future; for example, if the mild growth 
defect was fixed and both states grew at the same rate, the half- 
life of the ON state should be about 43 h. The sequestration- 
based switch should be advantageous for applications requir- 
ing scalability and predictable tuning. 

Materials and methods 
Modeling 

Bifurcation diagrams were computed from a simple model described in 
the Supporting information, using experimentally measured induction 
curves (shown in Supplementary Figures 9 and 10). The Supporting 
Information also provides an analytical treatment of the system, and 
the width of the hysteretic region as a function of the two key 
parameters is shown in Supplementary Figure 11. 



Cloning details 

The sigW and rsiW genes were cloned fromB. subtilis 168 (Bs) genomic 
template using primers o330/o331 for sigW and o332/o333 for rsiW 
(oligonucleotide sequences are in Supplementary Table 3). When we 
express the anti-sigma factor rsiW in Ec, we only use the N-terminal 87 
amino acids as this domain normally remains in the cytoplasm, does 
not contain residues that target RsiW for proteolytic degradation, and 
was previously found to be sufficient to sequester the sigma factor in 
Bs (Schobel et al, 2004; Zellmeier et al, 2006). For the SigW-dependent 
promoter, we cloned the promoter that drives the sigW/rsiW operon in 
Bs (using primers o328/o329) (Qiu and Hermann, 2001). We observed 
some leakage from this promoter in the absence of sigW expression, 
but this leakage was removed when we changed the UP element of the 
promoter to the Ec consensus (using oDC383/oDC384 and the Phusion 
Site-Directed Mutagenesis protocol). This modified promoter exhibits 
a 100-fold induction ratio and was used to control transcription of both 
the positive feedback loop and the YFP reporter. 

As explained in the Supporting information, increasing maximal 
sigma factor production rate from the feedback loop (QUOTE in the 
model) increases the width of the hysteretic region and increases the 
difference between ON and OFF states. We selected for high expression 
by constructing an RBS library in the 5'-UTR region between PsigWand 
sigW (using primers o362/o363) and screened for candidates that 
exhibited strong fluorescence when co-transformed with a reporter 
plasmid containing PsigW-GFP. We also created an RBS library for rsiW 
expression (using primers o364/o365) and selected for an RBS 
that produced enough anti-sigma factor to turn off fluorescence from a 
strain containing PsigW-GFP and sigW being expressed from the Ptet 
promoter. 

The PsigW-sigWand Pbad-rsiW transcriptional units were placed on 
a high copy plasmid (pMBl origin) with a strong transcriptional 
terminator (rrnB terminator) between them. The PsigW- Venus and 
Ptet-sigW transcriptional units were placed on a pi 5a plasmid, again 
with a strong transcriptional terminator between them. To facilitate 
debugging of the inducible promoters, RFP was added to the operon 
controlled by Pbad (using primers o404/o405 and o391/o406) and CFP 
was added to the operon controlled by Ptet (using o395/o396 and 
o390/o394) using SLIC (Li and Elledge, 2007). The Venus (YFP) and 
Cerulean (CFP) reporters were obtained from pZS2-123 (Cox et al, 
2010) . RFP is a codon-optimized mRFP obtained by Prof. J Christopher 
Anderson, UC Berkeley, from the Parts Registry (part E1010) and was 
modified at nucleotide positions 6-11 to fix an RBS issue. 

The sigE switch based on R. sphaeroides was cloned using o416/ 
o417 off pBS16 (Schilke and Donohue, 1995). The sigE/chrR system 
was previously shown to function in Ec (Newman et al, 2001). The 
reporter plasmid pDC311 was created by using o418/o419 off pDC288 
to swap the sigW-responsive promoter with the sigE-responsive PI 
promoter. Following the same procedure described above for the sigW- 
based switch, an RBS library was created between the PI -sigE 5'-UTR 
region using o420/o421 and the candidates were screened using 
pDC311. The chrR gene was also cloned from pBS16 using o422/o423 
and used to replace rsiW in pDC298 using SLIC with o424/o425. The 
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candidate from the RBS library (using o428/o429) and the Pbad-chrR 
(using o426/o427) operons were combined using SLIC to make 
pDC317. 

Plasmid maps are provided in Supplementary Figure 8 for the SigW- 
based bistable system and Supplementary Table 2 contains the 
plasmids and strains. Plasmids were transformed into BW27783, a 
strain designed to allow homogeneous expression from the pBad 
promoter (Khlebnikov et al, 2001). Oligos were ordered from IDT with 
standard desalting. 



Growth conditions 

Cells were grown at 37°C in M9 minimal medium supplemented with 
glycerol and casamino acids; media contains per liter: 11.28 g 5 x M9 
salt (Sigma M6030), 300 mg thiamine hydrochloride, 8ml 50% 
glycerol, 2g casamino acids (BD Bacto), 20 ml of 0.1 M magnesium 
sulfate, 200 ul of 0.5 M calcium chloride, filling to 11 with distilled 
water, and filtering before storage at 4°C. Antibiotics were added as 
appropriate with carbenicillin at 50ug/ml and chloramphenicol at 
25 ug/ml. To prepare master cultures, cells were grown overnight in 
0.1% arabinose (Sigma A3256) or the absence of arabinose to an 
OD600 of 0.1 and then stored in 10% glycerol aliquots of 50-100 ul at 
- 80°C (8-strip PCR tubes, T-3135-1, from Bioexpress were found to be 
the most convenient for storing aliquots). Master cultures for 
Figure 3B-D were grown in the absence of aTc and 0 or 0.1% 
arabinose. The experiments were performed by thawing a frozen 
aliquot at 37°C for a minute and then diluting the cells 16-million-fold 
into 2 ml of the m9 medium. Cells were grown overnight at 37°C at 
200 r.p.m. until an OD600 of ~ 0. 1 and then cells were concentrated via 
a microcentrifuge for microscopic analysis. For the tunable boundary 
experiments, cells were grown at the indicated concentrations of 
arabinose and aTc (Clontech 631310). For the set-reset latch experi- 
ments, cells were pulsed OFF by growing in 0. 1 % arabinose and in the 
absence of aTc. Cells were grown in the memory state in the absence of 
aTc and 6.3 * 10 ~ 5 % arabinose. Cells were pulsed ON by growing at the 
same arabinose concentration and ~ 10 ~ 3 uM aTc. 



Microscopy 

Cells were imaged on a Zeiss Axio Observer Dl using a 63 x plan- 
apochromat Ph3 oil-immersion objective. Cells were imaged with 
phase contrast and with the appropriate fluorescence filter set (38 for 
GFP/YFP, 45 for RFP, and 47 for CFP) with the following exposure 
times (150ms for RFP, 50ms for GFP/YFP, and 150ms for CFP). 
Fluorescence excitation was provided with an X-CITE 120XL. Images 
were captured with an AxioCam MRm using the provided Zeiss 
Axiovision software. 



Image analysis 

Images were analyzed using Matlab software, starting with the 
routines provided with CellTracer (Wang et al, 2010). Phase contrast 
images were first thresholded and then segmented using the 
CTiterativeSelectiveSegmentation function, which recursively sepa- 
rates the cells in an area of the image. Cells were then discarded if they 
failed checks on eccentricity or size (removes clumps of cells that 
could not be segmented as individuals). Sample segmentation results 
are presented in Supporting Figure 12. After segmenting the phase 
contrast images, the fluorescence images were first flat-field normal- 
ized. Then the fluorescence for each cell in the image was calculated. 
For each cell, the background fluorescence in a 20-by-20 pixel box 
surrounding the cell was calculated and then subtracted from the mean 
fluorescence in the area within the cell. Finally, dead cells or debris 
misclassified as cells were removed by excluding cells that failed to 
exhibit any CFP when expected (i.e., cells containing pDC304 and 
grown in aTc induction or cells containing pDC297 with CFP under 
constitutive expression) . For cells not expected to express CFP, cells 
were excluded if they failed to exhibit any RFP fluorescence. For most 
images, no more than 1 % of the cells after segmentation were excluded 
as dead. 



Other instruments 

Bulk fluorescence measurements were performed on a Tecan Infinite 
M1000 (RFP settings were 584/607 nm for excitation and emission 
with 5 nm bandwidths, 0 |is lag, 20 us integration, 50 reads, calculated 
z-position, and gain set to 170) or a Tecan Safire 2 (YFP measured with 
510/528 excitation and emissions wavelengths and 5 nm bandwidths, 
130 gain, 10 reads, and 20 (is integration). Optical density measure- 
ments were performed on an UVmini-1240 (Shimadzu). 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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